A Discrete-Event Analytic Technique for Surface Growth Problems 
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We introduce an approach for calculating non-universal properties of rough surfaces. The tech- 
nique uses concepts of distinct surface-configuration classes, defined by the surface growth rule. The 
key idea is a mapping between discrete events that take place on the interface and its elementary 
local-site configurations. We construct theoretical probability distributions of deposition events at 
saturation for surfaces generated by selected growth rules. These distributions are then used to 
compute measurable physical quantities. Despite the neglect of temporal correlations, our approxi- 
mate analytical results are in very good agreement with numerical simulations. This discrete-event 
analytic technique can be particularly useful when applied to quantification problems, which are 
known to not be suited to continuum methods. 
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I. INTRODUCTION 

Broad interest in surface growth problems and the in- 
terface motion stems from many applications that these 
problems have in practically all sciences. Some well- 
known examples are crystallization problems, epitaxial 
growth, wetting phenomena, electrophoretic deposition 
of polymer chains, and growth of bacterial colonies. The 
mainstream studies of physical mechanisms that lead 
to growth and roughening, and interface properties, fo- 
cus on the so-called universal properties, i.e., large-scale 
properties of growing surfaces as determined by univer- 
sal scaling exponents. A central point here is a coarse- 
grained description of the interface, which leads to an 
effective growth equation, typical for a universality class. 
It is safe to say that universal properties of time-evolving 
surfaces are well understood^. There are many sit- 
uations, however, where non-universal properties, i.e., 
those pertaining to the microscopic structure of the in- 
terface, are of importance in understanding the physical 
system. One case is the density of local minima of a 
virtual-time interface, evolving in parallel discrete-event 
simulations 2 - 3 - 4 . Another case is the movement of ran- 
dom walkers on a growing surface^. Another case is the 
density of maxima after solidification and processing of 
a surface, where this density is related to the frictional 
forced. Even for a one-dimensional interface, as in the 
growth of a two-dimensional crystal^, a general descrip- 
tion of the microscopic structure is challenging, both for 
analytics and for simulations, due to a variety of length 
and time scales. However, for a strongly non-equilibrium 
interface (e.g., in a cold system) a net picture of the in- 
terface gains and losses can be simplified. This leads to 
a variety of simulation models such as, e.g., ballistic de- 
position, Eden or solid-on-solid models. 

Here we introduce a discrete-event analytic technique 
that provides a means for calculating some non-universal 
properties of interfaces. The key idea is the construc- 
tion of a probability distribution for events that take 
place on the surface. This new approach gives a mean- 



field like approximation of averages that are determined 
in simulations. In Sec. [HJ we explain the main concept 
by examples of applications to selected non-equilibrium 
one-dimensional model interfaces. The method is sum- 
marized in Sec. IIIII 



II. THE METHOD 

We focus on a model in which a one dimensional (ID) 
surface grows by the deposition of local height increments 
r/k at L lattice sites, with periodic boundary condition 
h,L+i = hi (k enumerates the lattice sites, hk is the local 
height). The deposited rjk is a real positive number that 
can take on continuous values in the interval [0; oo). Our 
discussion is not specific to the probability distribution 
from which rjk is sampled. We consider a rough surface 
at saturation (Fig. [TJ. Such a surface can be generated, 
e.g., in the single-step solid-on-solid models^ or in ran- 
dom deposition with surface relaxation^ or by depositing 
random rjk at local surface minimum^. 

In the latter case, the growth is simulated by the 
nearest-neighbor rule: hk(t+l) = hk(t)+rjk(t), if /ifc(i) < 
mm{hk-i(t), hk+i(t)}-, otherwise, h k (t + 1) = h k (t) (t is 
the time index) . The above deposition rule produces sur- 
faces characterized by the universal roughness exponent 2 
a = 1/2. Our objective is to show that it is possible to 
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FIG. 1: A typical model surface in ID, L — 100. 
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FIG. 2: The four groups of elementary local surface con- 
figurations at the fc-th site. In accordance with the adopted 
growth rule, local A-configuration represents four types of 
discrete deposition events, and local configuration B (or C) 
represents two types of events. 



construct a theoretical probability distribution of depo- 
sition events on the surface generated by the above rule. 
This distribution can then be used to compute measur- 
able quantities at saturation, e.g., the mean density of 
local minima or maxima. Also, it can serve as the depar- 
ture point for the construction of other distributions that 
describe more complex growth processes. Our deriva- 
tion makes the following two simplifying assumptions. 
First, we neglect correlations between nearest- neighbor 
local slopes, which depend on the type of deposition, i.e., 
the distribution from which 77^ is sampled. Second, we 
neglect temporal correlations among the groups of the 
surface configuration classes, as explained later. Because 
of the above simplifications, our theoretical results for 
the averages are a mean-field like approximation to the 
averages measured in simulations. 

The adopted growth rule provides the principle for the 
classification of local lattice-site configurations of the in- 
terface. There are only four groups of these configura- 
tions (Fig. [21 • Each group corresponds to one of the 
four mutually exclusive discrete events A, B, C and D 
that take place at site k during the deposition attempt 
at t: "A" denotes an event when the deposition rule 
is satisfied from both sides, i.e., hf-(t) < hk-i(t) and 
hk(t) < hk+i(t); "B" denotes an event when the rule is 
not satisfied from the right; "C" denotes an event when 
the rule is not satisfied from the left; "D" denotes an 
event when the rule is not satisfied from either side. The 
above events can be mapped in unique way onto elemen- 
tary local site configurations A, B, C and D, presented 
in Fig. For a closed chain of sites, the set of events 
when all sites are simultaneously in the A-configuration 
is of measure zero, and not all L sites can have the same 
elementary configuration. Therefore, in the set of L sites 
there must be at least one site with configuration A. We 
assign the index k = 1 to one of the sites in the A- 
configuration and enumerate the other sites accordingly, 
progressing to the right. Its right neighbor (site k = 2) 
can be only either in configuration C or D. Similarly, its 
left neighbor (site k — L) can be only either in B or D. 



If site k — 2 is in C then site k — 3 can be only either in 
C or D. If site k — 2 is in D then site k = 3 can be only 
either in B or A. If site k — 3 is in B then its right neigh- 
bor can be either in B or A, etc. Starting from the site 
k = I and progressing to the right towards k — L, ap- 
plying these elementary neighbor rules, we can construct 
all possible configuration equivalency classes of the en- 
tire surface generated by the adopted growth rule. These 
can be categorized into groups (called p-groups), based 
on the number p of the deposition events on the surface 
at t, i.e., here, the number of local minima in the surface 
configuration (coded by "A") at t. The probability dis- 
tribution f(p; L) of the deposition events is obtained as 
the quotient of the multiplicity M(p) of the p-group and 
the total number M of configuration classes. (We omit 
t in the notation since the analysis concerns steady-state 
evolution.) The mean density (u(L)) of local minima is 
given by the generally valid formula: 

(«(L)) = 5>(p)/(p;L), (1) 
p 

where the summation is over all p-groups of the admissi- 
ble surface-configuration classes, u(p) is the density char- 
acteristic for each group, and f(p;L) is the frequency of 
the occurrence of thep-group at saturation, i.e., the prob- 
ability of generating a surface with p local minima. 
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FIG. 3: Binary tree for the construction of all possible con- 
figuration classes for L=5. 

For example, the binary tree for the construction of 
possible surface-configuration classes for L = 5 is shown 
in Fig. |3 Looking along its branches, starting from 
the leading A at the fixed k = 1 position, one identi- 
fies M — 8 configuration classes of the entire surface: 
(I) ACCCD; (2) ACCDB; (3) ACDBB; (4) ACDAD; 
(5) ADBBB; (6) ADBAD; (7) ADACD; and (8) ADADB. 
Each surface configuration represents a class of infinitely 
many topologically equivalent deformations (Fig.Q since 
the deposited random height increment is a real posi- 
tive number. There are only two p-groups. In the first, 
there are four classes with one letter A: M(l) = 4, 
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/(1;5) = 1/2, and u(l) = 1/5. In the second, there are 
four classes with two letters A: M(2) = 4, /(2; 5) = 1/2, 
and u{2) = 2/5. By Eq. (jTJ, (tt(5)) = 3/10. 






(6) 



FIG. 4: The graphs of possible surface-configuration classes, 
corresponding to the configurations read along the branches 
from Fig. (see text). Each graph represents a class of in- 
finitely many topologically equivalent deformations. 

To hnd /(p; L), one can exploit the recurrent structure 
of the binary tree for general L in counting classes that 
contain the A-configuration at exactly p sitesi Counting 
gives: M = 2 L ~ 2 , M(p) = (L - l)!/((2p - 1)\{L - 2p)\), 
and p = l,2, [L/2] ([L/2] = L/2 for even L; [L/2] = 
(L-l)/2 for odd L). Thus, 



f(p;L) 



i 



2 i-2 



L - 
2p 



(2) 



In deriving f(p; L) the underlying assumption is that 
at saturation any class of the entire surface configura- 
tions is generated with probability 1/M. The theoretical 
f(p; 50) is compared with the simulated distributions in 
Fig. [SJ We considered three deposition types: the Pois- 
son type, where r]k was sampled from the Poisson distri- 
bution; Gaussian type, where rjk was sampled from the 
Gaussian; and, the uniform deposition with rjk being a 
uniform random deviate. In all cases, Eq. (J2J) is a very 
close approximation to time-averaged distributions mea- 
sured in simulations. 

Having Eq. J2J, all moments of f(p;L) can be com- 
puted exactly^, e.g., for L > 4 its variance is 



a 2 (L)=J2(m-{p(L))) 2 f(m;L) = 



L-l 
16 



For L > 3 the mean number of local minima is 

L + l 



(p(L))=J2^f(m;L) = 



(3) 



(4) 



The mean density of local minima, presented in Fig. |HJ 
is simply obtained from Eq. (0J as (u(L)) = (p(L))/L. 

Suppose, the adopted growth rule is slightly modified 
so the deposition at a local minimum (site A) happens 
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FIG. 5: Probability distributions for L = 50: theoretical 
(histogram), and from simulations with Poisson (circles) and 
Gaussian (diamonds) depositions. Error bars give one stan- 
dard deviation from the mean time sequence at saturation. 
The frequencies were measured in 2048 independent trials. 
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FIG. 6: The steady-state mean density (u(L)) of local min- 
ima vs lattice size L. The curve is the analytical result 
(L + 1)/4L. The symbols present simulations with the Poisson 
(circles), Gaussian (diamonds) and uniform (squares) deposi- 
tions. The error bars are smaller than the symbol size. 



with probability qA- Then, the new probability distribu- 
tion Fiji; L, qa) is the product of the probability f{nA\ L) 
of drawing a surface that has exactly ua local minima, 
given by Eq. |J2J, and the probability B(n; ua, qA) that n 
depositions occur on this surface. The latter is given by 
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binomial distribution 



( n A 



B{n;n A ,q A )=^)q n A q n A A - v 



(5) 



where first factor gives the number of ways in which 
n depositions can be distributed among n A sites, and 
q A = 1 — q A is the complementary probability of the no- 
deposition event at site A. In the old rule: B(n; n A , q A = 



1) 



so, trivially, F(n\ L; q A ) = f(n A ;L). 



Consider next a more complicated rule, where at 
the i-th deposition attempt, the deposition of %(t) 
is not exclusively confined to local minimum but may 
also happen at other sites B, C or D with the cor- 
responding probabilities q B , qc an d qD, respectively. 
We assume now q A = 1. The probability F(n;L) = 
F(n;L;n B ,nc,n D \q B ,qc,qD) of having n = n A +riB + 
nc + n D deposition sites in the interface is the product: 

F(n;L) = f(n A ; L)B{n B ; N B ,q B ) 

x B(n c ;N c ,qc)B(n D ;N D ,q D ). (6) 

In Eq. ©, f(n A ;L) is the probability of generating a 
surface with n A local minima, each of which is the de- 
position site; the binomial B(n B ] N B ,q B ) is the prob- 
ability of having n B depositions at N B sites in the lo- 
cal B-configuration on the generated surface; similarly, 
B(nc', Nc, qc) an d B{njj\ Njy 7 qo) are probabilities that 
nc and n B depositions occur at Nc and Njj sites in the 
local C- and D-configurations, respectively. For any sur- 
face, N A + N B + Nc + Nd = L. For surfaces growing 
on a closed ring of L sites, the number N B of local max- 
ima matches the number N A of local minima. This gives 
N B + Nc = L — 2n A . By Eq. (Q, the mean number 
(n(L)) of deposition sites in the interface is 

[#J N B N c n A 

( n ( L )) = E E E E nF(n;L). (7) 

ua — 1 riB—0 nc—0 ric—0 

The mean density of these sites is (U(L)) = (n(L))/L. 

In Eq. 0, the summation over no is easily performed: 
J2no=o( n A + n B + n c + n D )B(n D ;N D ,q D ) = n A (l + 
qo) + n B + n C - This gives 



(n(L)) = f( n ^ L ) i n A (1 + qo) + S BWC ) , (8) 

where S B yc is the mean number of deposition events in 
the set of sites other than A or D. The computation of 
S B vCi though elementary, is a bit tricky since one must 
distinguish between depositions in the set of B-sites and 
depositions in the set of C-sites to avoid double counting: 

N B N c 

SbvC = E B ( n s;N B ,q B ) 2J B(n c ;N c ,qc) 

riB—O nc—0 

x (n B + n c ) = (L- 2n A )q B q c . 
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FIG. 7: The steady-state mean density (U(L)) of deposition 
sites vs lattice size L when q depends on integer N > 2: ana- 
lytical curves of Eq. 1101 and the data measured in simulations 
with the Poisson deposition (symbols). 



Substitution to Eq. (JSJ) gives for L > 3: 

( n ( L )) = 1 L + 1 . (9) 

where Q = 1 - Q, Q = q B + q c - q B q c - 

One application of Eq. © is the case when q B ~q c = 
q/2 and q D = q, where the probability q — \/2jN is 
parametrized by an integer N > 2. Here, the mean den- 
sity of lattice sites where the deposition takes place is 



<fJ{L)) = [l-\ 



1 qL-l 
4 L 



(10) 



Predictions of Eq. (JTJ3J, compared to simulations with 
the Poisson deposition in Fig. [7| are consistent with the 
fact that the approximation f(n A ;L) gives slight over- 
estimates of the density of local minima at the high-L 
end (Fig.|SJ. Nevertheless, the analytical (U(L)) reflects 
well the overall tendencies in the data. As q — > 0, (U{L)) 
approaches 1. In this limit-case each attempted deposi- 
tion event at each lattice site becomes certain. Such a 
case corresponds to the random deposition growth rule. 
The corresponding interface is entirely uncorrelated and 
never saturates. This is the asymptotic limit to which 
the curves in Fig. converge when N — > oo. 



III. OUTLOOK AND CONCLUSION 

In summary, we have introduced a new approach for 
the study of non-universal properties of surfaces. We 
showed the explicit connection between the interface 
morphology and the event statistics on the interface at 
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saturation. The key concept here is the ability to build 
distinct equivalency classes of the entire surface config- 
urations from its local-site configurations. In this con- 
struction, the growth rule defines the set of events that 
can be mapped in a unique way on the set of elementary 
local-site configurations. 

Small differences between simulation results and the 
approximate analytical results come mainly from ne- 
glecting temporal correlations among p-groups of surface- 
configuration classes. These correlations are intrinsically 
present in the computation of averages over time series in 
simulations but they are absent in our derivation. They 
depend on the probability distribution from which the 
deposited random height increments are sampled. The 
proper handling of temporal correlations is essential in 
the microscopic description of the growth phase (before 
saturation), which is a challenging problem still open. 
Some insight to this issue may be found in recent works 
of Rikvold and KolesikiS, who used the concept of equiv- 
alency classes in connecting the interface microstructure 
with its mobility for Ising and solid-on-solid models with 
various stochastic dynamics. 

Although the focus of the study presented here is on 
obtaining analytical results in the mean- field spirit, it 
is interesting to notice the sensitivity of simulation re- 
sults to the probability distribution used to select random 



height increments. Negligible variations between the de- 
rived analytical probability distribution and simulated 
distributions with Gaussian, Poissonian and uniform in- 
crements demonstrate that the theoretical distribution 
presents well the statistics of events on the surface. Tem- 
poral correlations between interfaces generated via ran- 
dom deposition at local minima deserve further analytical 
studies. 

In conclusion, we observe that our approach could be 
applied to a variety of other growth rules and models. 
The current applications to ID problems are illustrative 
examples. The main advantage of the approach is that it 
enables one to compute analytically quantities that can 
be only estimated qualitatively by continuum methods. 
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